Darcy Forchhiemer imposed exponential heat source-sink and activation energy with the effects of bioconvection over radially stretching disc

The Darcy–Forchheimer model is a commonly used and accurate method for simulating flow in porous media, proving beneficial for fluid separation, heat exchange, subsurface fluid transfer, filtration, and purification. The current study aims to describe heat and mass transfer in ternary nanofluid flow on a radially stretched sheet with activation energy. The velocity equation includes Darcy–Fochheimer porous media effects. The novelty of this study is enhanced by incorporating gyrotactic microorganisms which are versatile and in nanofluid can greatly improve the thermal conductivity and heat transfer properties of the base fluid, resulting in more efficient heat transfer systems. Furthermore, the governing PDEs are reduced to ODEs via appropriate similarity transformations. The influence of numerous parameters is expanded and physically depicted through the graphical illustration. As the Forchheimer number escalates, so do the medium's porosity and drag coefficient, resulting in more resistive forces and, as a result, lowering fluid velocity. It has been discovered that increasing the exponential heat source/sink causes convective flows that are deficient to transport heat away efficiently, resulting in a slower heat transfer rate. The concentration profile accumulates when the activation energy is large, resulting in a drop in the mass transfer rate. It is observed that the density of motile microorganisms increases with a rise in the Peclet number. Further, the results of the major engineering coefficients Skin-friction, Nusselt number, Sherwood number, and Microorganism density number are numerically examined and tabulated. Also, the numerical outcomes were found to be identical to the previous study.


Re
Local Now-a-days, the problems associated with boundary layer flow due to a stretching surface have gained prominence due to their practical utility in various industrial and technical workflows such as modern extrusion processes, extraction of copper wires, incremental sheet forming, and geothermal energy applications.The stretching rate and cooling rate have a major influence on the mechanical characteristics of the fluid that these systems demand as their output.As a result, the heat transfer properties of the stretching sheet have received considerable attention from researchers.Crane 1 was the first to work on a stretching sheet problem where he explored the boundary layer flow over the surface.Later on, a significant number of researchers reviewed the stretching sheet problems in the context of various fluid models.Cortell 2 addresses the numerical analysis for viscous flow and heat transfer across a non-linear stretching sheet.Reza-E-Rabbi 3 explored MHD nanofluid flow over stretching sheet, the author illustrated comparative between Casson and Maxwell fluid over different flow fields.Jeelani and Abbas 4 investigated suction and radiation effects over Maxwell hybrid nanofluid flow past a permeable inclined stretched sheet.MHD, radiation, and chemical reaction impacts over unsteady Casson nanofluid flow across a stretching sheet were investigated by Reza-E-Rabbi 5 .This article also explored the consequences of thermophoresis and Brownian motion over fluid flow.Considering the influence of radiation and magnetic field, Srinivasacharya and Kumar 6 studied Casson flow on a stretching surface.
The solution to the flow problem was computed using artificial neural networks, and it was discovered that this neural network method was accurate, and the effectiveness of the solution improved as the number of neurons in the neural network expanded.Yousuf Ali 7 examines the combined effects of Hall current and radiation on MHD nanofluid flow through a nonlinear stretching sheet.This study has found that when the effects of radiation and Hall current are combined, heat transfer is relatively higher than when only one of those effects is taken into account.With convective boundary conditions, Rafiqul Islam 8 examined the heat and mass transport phenomena of a Casson nanofluid flow under the influence of MHD, heat source/sink, and chemical reaction.
The existence of a heat source or sink can have a substantial impact on the fluid flow and heat transfer characteristics.Several studies have been conducted to study the impacts of heat source/sink on various fluid flow scenarios, including MHD flow, nanoparticle aggregation, divergent/convergent channels, and non-Newtonian fluid flow [9][10][11] .These studies show that the existence of heat sources or sinks could manipulate temperature distribution, velocity profiles, and other flow parameters.Abbas et al. 12 analyzed the implications of radiation and inclined MHD on nanofluid flow, taking into account exponential heat source/sink.Thermal radiation and chemical reactions are also investigated to determine the thermal and mass transfer properties of nanofluids.Khan et al. 13 explored the Darcy-Forchhiemer flow of an Eyring-Powell nanofluid exposed to an exponential heat source/sink and gyrotactic microorganisms.The Cattaneo-Christov theory equation is utilized to highlight the heat and mass transfer phenomenon.
Flow owing to the Darcy-Forchheimer medium is a significant aspect of industrial applications such as geothermal energy production, catalytic converters, oil recovery processes, and gas turbines.There is a vast literature available on this topic such as: in the context of Darcy-Forchheimer model, Colak et al. 14 investigated the Powell Eyring nanofluid over a stretching surface with bioconvection and artificial neural network.Mandal and Pal 15 found dual solutions for the convective-MHD hybrid nanofluid in the Darcy-Forchheimer porous medium on a decreasing surface.Furthermore, the flow and heat transfer phenomenon were thoroughly investigated using stability analysis and entropy generation.Using the MHD Darcy-Fochheimer model, Joshi et al. 16 investigated the effects of suction/injection and dissipation on hybrid nanofluid flow across a permeable stretched sheet.Moreover, for the blowing region, velocity profiles exhibit dual behavior as Forchhiemer and porosity parameters are increased.Paatanaik 17 examined the combined impacts of Brownian and thermophoresis processes on nanofluid flow behavior across stretched surfaces, which also takes MHD, heat sources, and chemical reactions into account.
Activation energy is also used to investigate the nature of reactants and the effect of catalysts on reactions.Because activation energy is critical in boosting the rate of chemical processes, it is helpful in a variety of realworld applications such as match igniting, fire suppression, enzyme action, and many others.Certain reactions have been known to proceed slowly or not at all, even in the absence of a catalyst.Thus, an insignificant quantity of energy, known as activation energy, is required to start a chemical reaction.Babu and Sathian 18 assessed the effects of activation energy and lower viscosity on water flow through CNTs.It was discovered that nonlinear fluctuation in viscosity increased the rate of fluid flow through carbon nanotubes.The bio-convective hybrid nanofluid flow across a riga plate in the presence of heat absorption and the activation energy was studied by A. Algehyne et al. 19 .However, many studies on activation energy have been carried out recently [20][21][22] .
Bioconvection is a phenomenon caused by the upward movement of microorganisms, which are denser than water, causing the topmost part of the liquid to become unstable and triggering the emergence of convection patterns.The movement of microorganisms in bioconvection can be purposefully regulated to produce certain outcomes, making it an effective way in biotechnological systems.Bioconvection has been established in studies to be essential for the development of medicine delivery systems, biofuel production, and biological polymer synthesis.Platt's reports 23 were the first to discuss the term bioconvection, where he discovered Bioconvection patterns in Cells of Free-Swimming Organisms.Lui et al. 24 investigated the numerical bio-convective evaluation for a rate-type nanofluid that was affected by unique slip characteristics and Nield thermal limitations.These recent studies [25][26][27] have explored Bioconvection for different fluids with porous, non-linear radiation effects, and MHD.The below table shows the novel work in comparison with surveyed literatures (Table 1).
High thermal conductivity, improved stability in base fluids, crystallinity, zero potential, and a huge surface area are the only benefits of ternary nanofluids.The most recent flow model has not yet been explored, as found by initial research.This study examines the Darcy-Forchheimer model with gyrotactic microorganisms of ternary nanofluid over a radially stretching sheet with activation energy.Radial stretching techniques have the potential to create unique copper wire products with specific geometries or performance requirements, such as those required for specialist electrical components or sophisticated manufacturing processes.This research helps in studying how usage of radially stretching sheets in the which makes copper wires suitable for electronic

Mathematical modeling
Consider Darcy-Forchheimer bio-convection flow of an incompressible nanofluid flow comprising nanoparticles (Cuo, AA7072, AA7075) with two base fluids ethylene glycol and water across a disk at z = 0 .The surface is elastic in nature and is stretching radially with velocity u w = ar , (a > 0) is constant as shown in Fig. 1b.Furthermore, exponential heat source/sink and activation energy are incorporated.T w and T ∞ , respectively represents the fluid's constant temperature at the wall and ambient temperature.The constant concentration and motile microorganisms on the sheet's wall surface are denoted by C w and N w , respectively, and C ∞ , N ∞ are the concentration and motile organisms at the far field.
The basic equations for the described problem can be expressed in the vector form as follows: (1) www.nature.com/scientificreports/Here, V is a function of (u, w) , and Ṽ = −bW c (Cw−C∞) ∇C .Under these essential stated assumptions, the PDEs form of the steady initial governing equations is expressed as per 28,31,32 .
with consistent boundary conditions where (u, w) are the velocity components along (r, z) axis.The term ν(µ/ρ) and α k/ ρCp denotes kinematic viscosity and thermal diffusivity, respectively.C b is drag coefficient, K is the permeability of the porous medium, Q * is exponential heat source, k r denotes chemical reaction parameter, E a is activation energy component and K is Boltzmann constant, and (T/T ∞ ) n e −E a /KT denotes modified Arrhenius function, b represents chemotaxis constant, D N represents Microorganisms Diffusion coefficient, and W c represents maximum cell swimming speed.Thermo-physical properties of base fluids (water and ethylene glycol) and nanoparticles are assumed to be unchanged with temperature, and are listed in Table 2.

Particles
with The variables included in the equations are as follows: 5 , . The physical quantities are given by, From the above equations the Skin friction, Nusselt number, Sherwood number, and motile density are obtained as follows:

Methodology
This section discusses the methodology used to deduce solutions as well as code validation.The equations are initially modeled as PDEs and then transformed into ODEs using similarity variables.The "bvp-4c function" (a built-in package in MATLAB) is used to deduce the solution of Eqs. ( 14)- (17), as well as the BCs (18).With the aid of the "three-stage Lobatto IIIA formula", the "bvp4c function" employs a finite difference scheme and fourth-order accuracy.To achieve accurate results, the suitable boundary layer thickness, point depicting the far field η ∞ , and initial guess must be selected based on the parameters applied.To determine the model's solution, we convert the revised equations into a first-order system by introducing additional parameters.By utilizing these new parameters, the Eqs.( 14)-( 17) are reduced by following substitutions (13)    www.nature.com/scientificreports/ The following MATLAB syntax is used for equations and BCs. and The numerical solutions are then derived by feeding Eqs. ( 24) through ( 27) into the bvp-4c solver.The "sol" component is the outcome of bvp4c."sol = bvp4c (@OdeBVP, @OdeBC, solinit)" provides the solver's syntax, which is made up of a couple of functions.The codes for Eqs. ( 24)-( 27) are encoded in the ''@OdeBVP'' function.The BCs are defined using the ''@OdeBC" function (28).The ''solinit" function is used to code the initial mesh points and the initial solution approximation at those places.The solver will then execute which requires the substitution of the parameters Fr = 0.2 , Pr = 6.3(water)/7.56Engine oil , Qs = 0.3 , m = 0.01 , η = 3 , Sc = 0.6 , Rc = 0.5 , δ = 0.1 , E = 0.1 , Pe = 0.5 , and Lb = 0.5 with characteristics listed in Table 1. and properties of ternary nanofluids.A graphical representation of the results achieved for each constraint by altering each parameter while keeping other parameters still and the results will be shown via graphs.Furthermore, the values of missing conditions are estimated to begin the process of finding the solution, and other parameters in Eqs. ( 24)-( 27) are set to get the desired result.The solution is accepted after iterating and satisfying the conditions in Eq. ( 28) asymptotically.Figure 2 depicts the flow chart for determining the solution.We have taken the mesh size equal to 100.The convergence of the solution depends on the η → ∞ values such as the values of the different parameters choosen by asymptotic convergence of the numerical outcome.www.nature.com/scientificreports/Furthermore, when the numerical outcomes were compared to those from previous studies (see Table 3), the results were quite matching.

Results and discussion
This section of the work addresses the thorough analysis and physical interpretation of the momentum, temperature, concentration, and motile density profiles of the (water and ethylene glycol)-based Cuo, AA7072, AA7075 ternary nanofluids.For various values of the relevant parameters encountered in the problem, the numerically generated results for the friction coefficient, heat and mass transfer coefficients, and motile density are plotted via graphs utilizing bvp4c solver.The current study's authors examined the model in two separate base fluid scenarios: water and engine oil.
The ***Table 4 shows the results computed of Engineering coefficients limiting certain variables.The range of the effective used for the computational purposes are as follows: 0 ≤ Fr ≤ 0.6 , 0.1 ≤ Qs ≤ 0.7 , 0.1 ≤ Rc ≤ 0.7 , 0.5 ≤ E ≤ 2.0 , 0.5 ≤ Pe ≤ 2.0 , 0.3 ≤ Lb ≤ 0.9 .Table 3 shows the results computed of f ′′ (0) , θ ′ (0) , χ ′ (0) , and γ ′ (0) with limiting certain variables.From the table it is clear that upsurge in the value of Fr increases surface drag coefficient f ′′ (0) and rate of heat transfer θ ′ (0) whereas drop in mass transfer χ ′ (0) and motile density γ ′ (0) is observed.The rise in the Qs has drop the thermal distribution by thickening of thermal boundary layer.Improve- ment in the values of Sc , δ and Rc will increase the χ ′ (0) and γ ′ (0) profiles.The decrease in χ ′ (0) is observed for rising values of E but the γ ′ (0) profile reduces due to decreased microorganism diffusivity in the fluid.The augmentation in the values of Pe and Lb increases γ ′ (0) enhancing swimming rate of motile microorganisms.Upsurge in the values of Pr has elevated f ′′ (0) and θ ′ (0) profiles but opposite trend is witnessed in case of χ ′ (0) and γ ′ (0) profiles.Table 3.Values of f ′′ (0) when M = 0 , Fr = 1 in the absence of A 1 and A 2 .

Sl. No
Previous study f ′′ (0) 36 Present study f ′′ (0) Table 4. Computational values of Engineering coefficients f ′′ (0) , θ ′ (0) , χ ′ (0) , and γ ′ (0). Figure 3a,b visualizes the influence of Forchheimer number (Fr) over velocity and temperature profile, respec- tively.A larger value of Fr number displays a drop in fluid's velocity whereas the opposite nature is seen in the temperature field for both cases i.e., water and engine oil.The explanation for this behavior is that the inertia of the porous media adds more resistance to the fluid flow, causing the fluid to move at a slower rate with a lower temperature.Due to their physical properties, in the velocity profile, the decline in Engine oil is more than in water but the temperature rise is more in water when compared to engine oil.
Figure 4 describes the enhancement in the temperature profile θ (η) under the impact of exponential heat source/sink parameter (Qs) .It is obvious from the figure that increase in Qs promotes temperature rises.The temperature rise is prompted by the transfer of heat to the system via the heat source's internal fluctuations.Due to varied thermal conductivity of two base fluids, the temperature elevation is higher in water rather than in engine oil.The implementation of exponential heat sources/sinks can optimize heat distribution and energy utilization in industrial processes like nuclear reactors, steel manufacturing, and many more.
Variations in χ(η) for different values of reaction rate parameter is observed in Fig. 5.The main reason for this is that when the Rc grows, the number of solute molecules undergoing chemical reaction gets bigger, result- ing in a drop in the concentration field.As a result, a damaging chemical reaction dramatically diminishes the solutal boundary layer thickness.Due to the variations in their chemical characteristics, chemical reactions in engine oil tend to proceed more slowly than in water.Hence the reduction in concentration is seen more in water than in engine oil.their chemical properties, elevations in concentration are witnessed more efficiently in engine oil than in water.In thermal research, activation energy is used to determine how much energy is required for molecules to shift from one state to another.This information helps to improve temperature-dependent reactions and thermal management systems.
The effect of Peclet number (Pe) against motile density profile γ (η) is portrayed in Fig. 7.The figure indicates that fostering Pe results in microbe density condensing.Physically, a faster rate of advective movement results in a greater Peclet number, which rapidly increases the flux of microorganisms.Considering the Peclet number is inversely related to cell swimming speed, boosting the Peclet number can result in a drop in the density profile of motile microbes.
Figure 8 exhibits the impact of Lewis number (Lb) on γ (η) .Escalating the values of Lb tends to diminish the γ (η) .From a physical point of view in a bioconvection flow scenario, the Lb regulates the behavior and diffusivity of microorganisms.A rise in the bioconvection Lewis number causes microbe diffusivity to go down resulting in a reduction in their profile.Bioconvection happens in both water and engine oil, but its factors and implications differ, therefore water has a greater decline in motile density than engine oil.
The Darcy-Forchheimer model is critical for optimizing thermal management systems.For higher values of the Forchhiemer number, the skin-friction coefficient increases, influencing the efficiency of the heat  www.nature.com/scientificreports/transmission process.Figure 9 provides useful insights into the effect of Fr on skin friction Cf for two different base fluids water and engine oil.From a physical perspective, the inertia coefficient is directly proportional to the porosity of the medium and the drag coefficient, as Cb grows, so does the porosity of the medium and the drag coefficient.As a result, the Cf is increased, resulting in a decreased velocity.Because of their physical qualities, engine oil has a greater increase in the Cf than water.
In some industrial processes, particularly those that require precise control over heat transfer rates, where excessive heat transfer is undesirable or to avoid overheating, minimizing the Nusselt number can help improve energy efficiency by eliminating unnecessary heat loss.Figure 10 provides considerable details on the contrast between the exponential heat source/sink Qs and Nu profile in the context of a ternary nanofluid flow, where the base fluids are engine oil and water.The Nu is estimated to decrease as the Qs goes up.The Nu drops because an increase in the aforementioned parameter induces a rise in the thermal boundary layer structure, resulting in less heat transfer rate.From the figure, it is very clear that the drop is significant in engine oil.
The impact of activation energy (E) on Sh is demonstrated in Fig. 11.It is clear from the figure that Sh curves rise for augmented values of E .From the physical point of view, the parameter E can influence the concentration distribution and mass transfer rate, which are able to alter the Sherwood number.A lower mass transfer rate, and hence a lower Sherwood number, could occur from a larger E values.The impact of E on concentration in engine oil decreases faster than in water due to differences in their chemical properties.The drop in the Sherwood number as activation energy grows reveals that convective mass transfer becomes more dominant than diffusion, resulting in modifications to the system's mass transfer characteristics.

Concluding remarks
Current research incorporates a novel study of TNF across a Darcy-Forchheimer medium over a radially stretching sheet in the presence of gyrotactic microorganisms.The flow is also accomplished by exponential heat source/ sink and activation energy.This study offers a comparative heat and mass transfer analysis between two base fluids engine oil and water which provides valuable data that engine oil excels in lubrication and heat dissipation for machinery whereas water is often used in cooling systems and in steam generations.The bvp4c approach is used for tackling the reduced ODEs.The major outcomes of the current problem are: (a) The fluid velocity reduces with increasing Fr values but opposite trend is witnessed in case of temperature profile.We were able to successfully interpret Darcy Forchiemer's medium effect over ternary nanofluid flow with exponential heat source/sink and activation energy across a radially expanding sheet using computational methods.This research could be developed to incorporate other generalized non-Newtonian fluid models, CNTs, Hall effects, and stagnation point flow over various geometries.

Figure 2 .
Figure 2. Flow chart of the problem.

Figure 6 Figure 3 .
Figure3a,b visualizes the influence of Forchheimer number (Fr) over velocity and temperature profile, respec- tively.A larger value of Fr number displays a drop in fluid's velocity whereas the opposite nature is seen in the temperature field for both cases i.e., water and engine oil.The explanation for this behavior is that the inertia of the porous media adds more resistance to the fluid flow, causing the fluid to move at a slower rate with a lower temperature.Due to their physical properties, in the velocity profile, the decline in Engine oil is more than in water but the temperature rise is more in water when compared to engine oil.Figure4describes the enhancement in the temperature profile θ (η) under the impact of exponential heat source/sink parameter (Qs) .It is obvious from the figure that increase in Qs promotes temperature rises.The temperature rise is prompted by the transfer of heat to the system via the heat source's internal fluctuations.Due to varied thermal conductivity of two base fluids, the temperature elevation is higher in water rather than in engine oil.The implementation of exponential heat sources/sinks can optimize heat distribution and energy utilization in industrial processes like nuclear reactors, steel manufacturing, and many more.Variations in χ(η) for different values of reaction rate parameter is observed in Fig.5.The main reason for this is that when the Rc grows, the number of solute molecules undergoing chemical reaction gets bigger, result- ing in a drop in the concentration field.As a result, a damaging chemical reaction dramatically diminishes the solutal boundary layer thickness.Due to the variations in their chemical characteristics, chemical reactions in engine oil tend to proceed more slowly than in water.Hence the reduction in concentration is seen more in water than in engine oil.Figure6reveals the physical effect of activation energy E on χ(η) .Larger values of E augments the χ(η) .The Arrhenius function deteriorates as the activation energy value increases rapidly, resulting in a stimulation of the generative chemical reaction and an improvement in the concentration field.Because of the differences in

Figure 4 .
Figure 4. Upshot of Qs over temperature profile.

Figure 5 .
Figure 5. Upshot of Rc over concentration profile.

Figure 6 .
Figure 6.Upshot of E over concentration profile.

Figure 7 .
Figure 7. Consequence of Pe over motile density profile.

Figure 8 .
Figure 8. Consequence of Le over motile density profile.

Figure 11 .
Figure 11.Outcome of Sh over Sherwood number.

Table 1 .
Comparison of the current study with previously published papers.